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We have measured the Faraday rotation toward a large sample of polarized ra¬ 
dio sources behind the Large Magellanic Cloud (LMC), to determine the struc¬ 
ture of this galaxy’s magnetic field. The magnetic field of the LMC consists of 
a coherent axisymmetric spiral of field strength ~ 1 microgauss. Strong fluc¬ 
tuations in the magnetic field are also seen, on small (<0.5 parsecs) and large 
(~ 100 parsecs) scales. The significant bursts of recent star formation and su¬ 
pernova activity in the LMC argue against standard dynamo theory, adding 
to the growing evidence for rapid field amplification in galaxies. 

The Milky Way and many other spiral galaxies show well-organized, large-scale magnetic 
fields (dEHHl, the existence of which points to a powerful and ubiquitous process which or¬ 
ganizes random motions into coherent magnetized structures. The underlying mechanism is 
believed to be a dynamo, in which magnetic fields are slowly ordered and amplified due to the 
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interplay between turbulenee and differential rotation (|?l|5l). Magnetism in galaxies is usually 
mapped by observing the orientation of polarized optieal and radio emission from the galaxy 
itself (0), but these data have limited spatial resolution © and ean be diffieult to interpret (0. 
An alternative, direet determination of the geometry and strength of magnetie fields eomes from 
the Faraday rotation of baekground radio sourees dSl 19117^. an effeet in whieh birefringenee in 
an intervening magneto-ionized souree rotates the plane of linearly polarized radiation. Mea¬ 
surements of baekground rotation measures (RMs) are free from the diffieulties assoeiated with 
studying polarized emission produeed by the souree itself, whieh suffer from a eomplieated 
eombination of internal and external Faraday rotation, depolarization, and optieal extinetion. 

We have studied the magnetie field of the LMC, using 1.4-GHz polarization data reeorded 
as part of a hydrogen line survey (1771) earned out with the Australia Teleseope Compaet Ar¬ 
ray dH). Over a field of 130 square degrees, we have ealeulated RMs for 291 polarized baek¬ 
ground sourees. Of this sample, about 100 sourees lie direetly behind the LMC. We used 140 
measurements of sourees lying outside the LMC to subtraet a mean RM, presumably resulting 
from foreground Faraday rotation in the Milky Way. The resulting distribution of residual Fara¬ 
day rotation (Fig. 1) shows a strong exeess in RM aeross the extent of the LMC. The implied 
magnetie field demonstrates spatial eoherenee: the eastern half of the galaxy shows predomi¬ 
nantly positive RMs, while in the west the RMs are mainly negative. 

We have eonverted position on the sky to a loeation within the LMC disk for eaeh RM 
measurement, assuming that the galaxy is inelined to the plane of the sky at an angle i = 35°, 
with its line of nodes at a position angle 0 = 123° (measured north through east) (ITib . The 
resulting dependenee of RM against position angle within the disk of the LMC (Fig. 2) shows a 
systematie variation, with the maximum mean RM oeeurring near the line of nodes. These data 
ean be well fit by a eosinusoid of amplitude RMq = 53 ± 3 rad m“^, demonstrating that the 
eoherent eomponent of the LMC’s magnetie field has an axisymmetrie spiral geometry (I7?b . as 
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is seen in other galaxies (□I^IT^b . The phase of this eosinusoid eorresponds to the piteh angle, p, 
of the spiral field dH, but in this ease we ean only infer a weak eonstraint, \p\ 20°, espeeially 

given that the uneertainty on 0 is > 10° (I7?b . 

In addition to the eoherent field, a strueture funetion analysis indieates random fluetuations 
in RM, with a standard deviation Srm = 81 rad m“^, and oeeurring on a oharaeteristie angular 
seale of ~ 0.1°, or L 90 parsees at the distanee to the LMC of 50 kpe. This may represent 
the evolved supernova remnants (SNRs) and wind bubbles whose interloeking shells dominate 
the morphology of ionized gas in the LMC on this seale 

To estimate the relative strength of the ordered and random field eomponents, we assumed 
that ionized gas in the LMC eonsists of a disk of projeeted thiekness D, in whieh eells of linear 
size L eontain elumps of ionized gas of filling faetor / and density Ue (17^ . In eaeh eell the 
magnetie field is eomprised of a uniform eomponent of strength Bq plus a randomly oriented 
eomponent of strength Br. If Bji is uneorrelated with fluetuations in Ue, it ean be shown that: 
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The oeeupation length, fD, of ionized gas is DMq/EMq ~ 530 pe, where DMq = / n^dl 
100 em“^ pe and EMq = f nidi ~ 19 pe em“® are the average dispersion measure (DM) and 
median extinetion-eorreeted emission measure (EM) integrated through the EMC, respeetively 
( 123 ). Equation dB then implies Br/Bq = 3.6: the random field dominates the ordered field, as 
seen in many other galaxies d^lITTIl. 

The aetual values of Bq and Br ean also be estimated. If Br and fluetuations in Ue are 
uneorrelated, then the strength of the ordered eomponent of the EMC’s magnetie field is 
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where K = 0.81 rad m“^ pe“^ em^ pG~^. The strength of the random field is then Br = 
3.6i?o ~ 4.1 pG, and the total magnetie field strength on large seales is Br = (Rq + ~ 
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4.3 yuG. We note that in seleeted regions where B and rie are correlated, as might result from 
compression in SNR shocks, the above approach overestimates Bq and underestimates Br, each 
by factors of ~ 2 (ITTIi . 

The polarized background sources are not randomly distributed across the LMC’s extent, 
but tend to avoid areas of bright Ha emission (Fig. 1). Specifically, background sources are 
increasingly depolarized as they propagate through regions of higher EM (Fig. 3). This suggests 
that we are observing beam depolarization, in which small-scale fluctuations in foreground RM 
produce interference between polarized rays along adjacent sightlines 0. Beam depolarization 
is only significant when the angular scale of RM fluctuations is smaller than the resolution of 
the data and the scale of intrinsic polarized structures. While the resolution here, ~ 40 arcsec, is 
comparatively large, extragalactic sources are typically of much smaller angular extent: for 1.4- 
GHz flux densities in the range 10-300 millijanskys as observed here, the median angular size 
is only 6 arcsec (dll, or 1.5 parsecs when projected against the EMC. The depolarization 
(Fig. 3) implies strong RM fluctuations on scales I 1.5 parsecs. In such a situation, beam 
depolarization reduces the intrinsic linearly polarized intensity, Pq, to a level (ITPIItIi: 

P = (3) 

where A is the observing wavelength and ctrm is the standard deviation in RMs across the 
source. To account for the dependence of P/Pq on EM (Fig. 3), we need to relate ctrm to the 
EM along a given sightline. If the magnetic field is uncorrelated with the ionized gas density, 
we expect that urm = k where /c is a constant. With this assumption, the fluctuating 

magnetic field on a scale I parsecs needed to produce the observed depolarization has a strength 
Br ~ yuG. In Figure 3, Equation © has been fit to the data for Pq ^ 0.104 and 

k ^ 1.8 rad m“^ pc“^/^ cm^. Assuming I < 0.5 parsecs, we find that Br > 5 yuG. We thus 
infer that there are significant RM and magnetic field fluctuations on sub-parcsec scales in the 
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ionized gas of the LMC. This phenomenon is also seen in our own Galaxy, and may trace the 
turbulent winds and H ll regions of individual stars (ESIEB- 

Most spiral galaxies are long-lived systems that exhibit significant rotational shear and that 
experience relatively constant star-formation rates over long periods of time. Coherent mag¬ 
netic fields in these galaxies are believed to be produced by a dynamo mechanism, in which 
small-scale turbulent magnetic fields are amplified and ordered by cyclonic motions and differ¬ 
ential rotation However, in galaxies dominated by sudden bursts of star formation and 

supernova activity, the dramatic injection of energy should disrupt the slow mono tonic increase 
of the large-scale field produced by a standard turbulent dynamo (|221E0. The LMC has ex¬ 
perienced several intense bursts of star formation over the past ~ 4 Gyr triggered by repeated 
close encounters with the Milky Way and with the Small Magellanic Cloud (I24II25I) . and yet still 
maintains a coherent spiral magnetic field. Combined with previous results demonstrating the 
presence of ordered magnetic fields in young galaxies for which the dynamo has had little time 
to operate (|2^ . and in irregular galaxies which lack significant amounts of rotation (23, there 
is now evidence that standard dynamo processes are ineffective in the LMC and these other 
galaxies. There are several viable alternatives to explain the coherent magnetic fields that we 
observe. Potentially most pertinent for the LMC is the cosmic-ray driven dynamo, in which re¬ 
cent supernova activity generates a significant population of relativistic particles. The buoyancy 
of these particles inflates magnetic loops out of the disk; adjacent loops reconnect, and then are 
amplified by differential rotation to generate a large-scale spiral field (l2<jilL29ll . This mechanism 
not only requires vigorous star formation, as has occurred recently for the LMC, but has a time 
scale for amplification of only ~ 0.2 Gyr (1^^ . and so can quickly generate large-scale mag¬ 
netic fields before they are dissipated by outflows and tidal interactions. This process can thus 
potentially account for the coherent fields seen in the LMC and other galaxies (ESI. 
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FIGURE CAPTIONS 


FIGURE I: Faraday rotation measures through the Large Magellanie Cloud. The image shows 
the distribution of emission measure toward the LMC in units of pe om~®, derived from the 
Southern H-Alpha Sky Survey Atlas (U7T) . The symbols show the position, sign and magnitude 
of the baseline-subtracted RM measurements dH): filled and open circles (both marked in 
green) correspond to positive and negative RMs, respectively, while asterisks (marked in purple) 
indicate RMs which are consistent with zero within their errors. The diameter of each circle is 
proportional to the magnitude of the RM, the largest positive and negative RMs being +247 ± 
13 rad m“^ and —215 ± 32 rad m“^, respectively. 

FIGURE 2: RM against position angle within the LMC. The six data points are a binned 
representation of the 93 RMs which lie within a radius of 3.5 degrees of the center of the 
EM distribution seen in Fig. 1 (i.e.. Right Ascension [J2000] 05^16™03®, Declination [J2000] 
—68°41'45"), plotted against deprojected position angle within the LMC, measured from the 
line of nodes. The uncertainty on each datum is the weighted standard error in the mean for 
RMs in that bin. The dashed line shows a cosinusoidal least-squares fit to the unbinned data, 
with an amplitude of +53 ± 3 rad m“^ and an offset from zero of +9 ± 2 rad m“^. The phase of 
the cosinusoid is only weakly constrained, falling between ±15°. The fit is not a strong function 
of the center adopted for the LMC. 

FIGURE 3: Polarized fraction of 81 background sources as a function of EM (12 of the 93 
sources shown in Fig. 2 have been excluded: six with EM > 100 pc cm“®, and six with an 
observed EM < 0 pc cm“® due to imperfect star subtraction). A Galactic foreground contribu¬ 
tion of 3 pc cm“® has been subtracted from each EM measurement. The uncertainty for each 


9 


binned data-point corresponds to the weighted standard error in the mean for each bin. The 
observed depolarization as a function of EM cannot be a result of source confusion or other 
observational selection effects, since sources with RMs were identified from an image of linear 
polarization, in which the weak signals from diffuse polarized emission show no correlation 
with Ha emission. It also cannot be due to excessive Faraday rotation across our observing 
band (bandwidth depolarization), since for the narrow frequency channels (8 MHz) used here, 
this effect would manifest itself only for |RM| > 4000 rad m“^, ~ 20 times larger than any 
RMs observed. The dashed line shows a least-squares fit of Equation Q to the unbinned data, 
assuming ctrm cc 
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SUPPORTING ON-LINE MATERIAL 


Materials and methods 

The RMs presented here were determined from a 1.4-GHz mosaic of 1300 pointings toward 
the LMC, made using the Australia Telescope Compact Array over the period 1994 October 
to 1996 February at a resolution of 40 arcsec (Kim et al, 1998, The Astrophysical Journal, 
v503, p674). Full polarization was recorded in this survey in 14 successive channels of band¬ 
width 8 MHz each, with center frequencies between 1328 and 1432 MHz. Antenna gains were 
calibrated using observations every ~ 30 min of one of PKS B0407-68 or PKS B0454-810, 
while absolute flux densities and polarization leakages were determined using observations of 
PKS B1934-638. For each frequency channel, images were formed in Stokes Q and U at a 
pixel size of 13 x 13 arcsec^, were deconvolved using a maximum entropy algorithm, corrected 
for primary beam attenuation, and then were combined and debiased to form an image of linear 
polarization, C = Using a false-discovery rate algorithm (Hopkins et al., 2002, 

The Astronomical Journal, vl23, pl086) on an average of the 14 images of C, 324 polarized 
sources were identified which were unresolved, did not correspond to a catalogued pulsar or 
supernova remnant within the LMC, and had a fractional linear polarization between 0.3 and 50 
per cent. For each such source, values of Stokes Q and U were then extracted from the individ¬ 
ual channel maps for nine pixels surrounding the source peak. For each pixel, the RM and its 
uncertainty were determined from the 14 pairs of Q and U values using the algorithm of Brown 
et al, 2003, The Astrophysical Journal Supplemental Series, vl45, p213. An individual pixel’s 
RM measurement was accepted only if the pixel had a debiased signal-to-noise in polarization 
of > 1 in at least eight frequency channels, and if the quality of the RM fit was better than 
90%. A source’s overall RM was considered valid if at least five of nine pixels had acceptable 
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RM fits, and if the dispersion in RM between pixels was less than twiee the average error in 
eaeh pixel’s RM. Of 324 polarized sourees, 291 passed all these tests, with typical errors in 
RM for each source of ±20 rad m“^. By considering 140 RMs at a radius of > 4.5 degrees 
from the LMC’s center, we found that a positive offset in RM was present due to foreground 
and background contributions; a mean baseline was therefore subtracted from the data in four 
separate quadrants: ±33.5 ±3.3 rad m~^ in the southeast, ±22.5 ±3.8 rad m“^ in the southwest, 
±49.0 ± 4.5 rad m“^ in the northwest, and ±31.3 ± 3.4 rad m~^ in the northeast. 

To estimate the total dispersion measure through the LMC, we used the fact that disper¬ 
sion of the pulsed signals from the five known radio pulsars in the LMC implies a mean total 
dispersion measure DM ~ 100 cm“^ pc (Crawford et ah, 2001, The Astwphysical Journal, 
v553, p367). The foreground contribution to the DM from the Milky Way is ~ 50 cm“^ pc 
(Cordes & Lazio, preprint, http://arxiv.org/abs/astro-ph/0207156). If the mean DM of LMC 
pulsars corresponds to a location half way through the LMC disk, then DMq ~ 100 cm“^ pc. 

Emission measures have been derived from the Southern H-alpha Sky Survey Atlas (Gaus- 
tad et ah, 2001, Publications of the Astronomical Society of the Pacific, vll3, pl326), by cor¬ 
recting smoothed and star-subtracted Ha emission in this region for both foreground extinction 
in the Milky Way and internal extinction in the LMC. In both cases, extinction corrections were 
derived from Hi column density maps (Staveley-Smith et ah, 2003, Monthly Notices of the 
Royal Astronomical Society, v339, p87; Kim et ah, 2003, The Astrophysical Journal Supple¬ 
mental Series, vl48, 473) using typical gas-to-dust ratios for our Galaxy and for the LMC (Pei, 
1992, The Astrophysical Journal, v395, pi30), and assuming that extinction within the LMC is 
caused by only half of the LMC H l column. Ha intensity was then converted into EM units 
assuming an electron temperature of 8000 K. 
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